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We propose a class of ground states for doped Mott insulators in the electron second quantization 
representation. They are obtained from a bosonic resonating valence bond (RVB) theory of the t — J 
model. At half filling, the ground state describes spin correlations of the S=l/2 Heisenberg model 
very accurately. Its spin degrees of freedom are characterized by RVB pairing of spins, the size of 
which decreases continuously as holes are doped into the system. Charge degrees of freedom emerge 
$H ■ upon doping and are described by twisted holes in the RVB background. We show that the twisted 

holes exhibit an off diagonal long range order (ODLRO) in the pseudogap ground state, which has 
a finite pairing amplitude, but is short of phase coherence. Unpaired spins in such a pseudogap 
ground state behave as free vortices, preventing superconducting phase coherence. The existence 
of nodal quasiparticles is also ensured by such a hidden ODLRO in the ground state, which is non 
Fermi liquid-like in the absence of superconducting phase coherence. Two distinct types of spin 
excitations can also be constructed. The superconducting instability of the pseudogap ground state 
is discussed and a d-wave superconducting ground state is obtained. This class of pseudogap and 
superconducting ground states unifies antiferromagnetism, pseudogap, superconductivity, and Mott 
physics into a new state of matter. 



> . I. INTRODUCTION 

An important question posed by the study of high temperature superconductors is whether the superconducting 
^vq | ground state of doped Mott insulators can be described by a BCS-like wave function. A d-wave BCS wave function 
and/or its derivatives have been widely used in the phenomenology of high temperature superconductors, partly 
| because electron pairing, as evidenced by 2e flux quantization and similar experimental results, is observed in the 
cuprate superconductors. However, owing to the effect of strong on-site Coulomb repulsion between the electrons, 
superconductivity in these systems may be quite different. Whereas in a conventional BCS superconductor, screening 
and retardation effects serve to minimize the role of Coulomb interaction, the role of the latter becomes crucial in 
a doped Mott insulator where the charge degrees of freedom are partially frozen by the severity of the Coulomb 
interaction. Strong arguments have been presented in the literature that a high temperature superconductor evolves 
from a Mott insulator doped by holes, and consequently the ground state can be very different from a conventional 
BCS statek 

A simple step towards incorporating Mott physics into a BCS description was suggested many years agoik This 
involves a projection (out of the Hilbert space) of doubly occupied sites in a <i-wave BCS stated. In this approach, 
• i-H , the charge degrees of freedom on a lattice site i, are partially frozen by a Gutzwiller operator (1 — n^n^). At half 
filling, only spin degrees of freedom survive, and the wave function describes a Mott insulating state. This class of 
wave functions has achieved some success in the phenomenology of the cuprate superconductors^^. A recent paper 
argues the case further—. However, the treatment of low energy spin degrees of freedom in this description is not 
satisfactory. Experimentally, strong antiferromagnetic (AF) correlations have been observed for the Cu spins, but 
projected d-wave BCS states cannot account correctly for such correlations. Since we opine that spin correlations are 
intimately related to superconductivity in a doped Mott insulator, we are motivated to seek an alternate description 
in which the spin degrees of freedom can be treated more accurately and systematically. 

In this paper, we propose a new class of ground states based on a microscopic treatment of the t—J model. These 
ground state wave functions are different from the projected BCS wave functions, and possess the following properties: 
(i) the spin degrees of freedom can fully restore the AF long range order (AFLRO) as well as low-lying spin wave 
excitations at half filling, and evolve into a spin liquid with strong short-range AF correlations at low doping; (ii) 
the no double occupancy constraint is always satisfied for any doping, and is incorporated self-consistently instead 
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of being enforced by a "brute-force" projection; (iii) there is an inherent d-wave superconducting instability at finite 
doping. Thus, these wave functions offer a unified description of antiferromagnetism and superconductivity in a doped 
Mott insulator. 

To motivate the form of the ground state wave function, let us first consider the ground state of the undoped Mott 
insulator on a square lattice. It is well known that the dynamics of the low energy degrees of freedom in this case 
are described by the S — 1/2 Heisenberg Hamiltonian. While the exact ground state of this model is not known, a 
comparison with numerical results shows that the best variational state, |\&o)j is given by 8 

l*o)= E W(i 1 -j 1 )...W(i n -j n )(i 1 j 1 )...(i n j n ) . (1) 

The ground state QJ is defined to be the bosonic RVB state of the insulatoA Here, (ij) stands for a singlet spin 
pairing between opposite sublattice sites i and j, and W(ij), the positive weight factor associated with it. The 
motivation for such a construction is to ensure that the ground state satisfies the so-called Marshall sign ruleS. 
Marshall showed that the ground state of the 5=1/2 Heisenberg Hamiltonian on a square lattice satisfies the 
condition, sgn [^(^li &2, ■■■, (T n)] — (— l) p ( CTl '°' 2 ' - ' <Tjv ). Here, P(cr 1 ,cr 2 , ...,0^) denotes, say, the number of down spins 
on sublattice A for the spin configuration |<7i,<72, ...,cjv). Then, the matrix elements of the Heisenberg Hamiltonian 
are negative and the ground state energy assumes the form of a sum over negative termsAS. Since each singlet bond 
(ij) satisfies the Marshall sign, it follows that the wave function of QJ, which is given by ^ Yl Tltij) ( — l) l W(i — j) 
[here i(j) stands for up (down) spin sites and h(i — j) = Q if ij S the same sublattice sites], always satisfies this 
criterion. 

The ground state (IJ describes both short-range and long-range AF correlations very accurately as the best varia- 
tional wave function for the insulator. It is then reasonable to expect that any putative RVB state, away from half 
filling, evolves from (IJ. With this in mind, we ask how the wave function QJ is modified in the presence of holes. As 
discussed above, singlet pairing is characterized in the wave function (IJ by the Marshall sign. It can be shown that 
the Marshall sign would be still obeyed if doped holes remain static. Then, from the form of (IJ, one may easily see 
that the nearest-neighbor (nn) hopping of holes, which displaces spins as backflow, generally leads to disordering of 
the Marshall sign (— l) p . Consequently, a hole moving along a closed path picks up a generalized Berry's phase, which 
was first identified by one of us and co-authors 11 . This effect, described (and called) as a phase string is illustrated 
in Fig. 1(a). It is a singular effect and cannot be repaired by low energy (transverse) spin excitations^. 



o I 
















/ 


+ 


+ / 


- 



Phase string effect 
(a) 



4~ 

s y •*« 



^ i 

■: i 

■■■ / + 

... f 

— rf'S 

■-.J -r 



spin-charge entanglement 
(b) 



FIG. 1: (a) The ground state picks up a sequence of signs, ( — 1) x (+1) x (+1) x known as a phase string, during an 
operation as a hole moves through a closed path and back to the origin with the displaced spin configurations being restored. 
This irreparable phase string effect has been shown to be a general consequence of the t — J model—, (b) The + or — sign in a 
phase string is determined by the exchange of the hole with an f or J, spin on the path, which can be viewed as an entanglement 
between the hole and a spin, illustrated by two entangled loops as both the hole and a displaced spin should return to their 
origins after the operation. 

It follows that a bare hole created in the RVB background QJ cannot possibly describe a low-energy state, and must 
be dressed by the nonlocal phase string effect. Since the phase string effect describes a mutual entanglement between 
spin and charge degrees of freedom illustrated in Fig. 1(b), the RVB background must also respond nonlocally. One 
may thus surmise that low-energy charge degrees of freedom are created by Q ff = Ci CT e~ je ", instead of by the bare Ci ai 
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where e~ l&a is a nonlocal operator describing the phase string effect. We call a hole that is created by Cj CT a twisted 
hole, and as will be shown in this paper, such a twisted hole is a bosonic object whose kinetic energy is optimized in 
the (modified) RVB background. We shall obtain a ground state of the following form 

(\ N l / \ N i 

5> T (0q t J (l>lW2uj IKVB), (2) 

where Z a (l) describes the wave function of the twisted hole q ct , and + = Nh is the total number of holes. The 
modified RVB vacuum |RVB) has the same structure as l^o) at half-filling, except that the RVB pairing amplitudes, 
Wij's, become complex and change continuously with doping, which can be determined self-consistently as we shall 
show. 

The ground state J2J is one of the central results of this paper. It generalizes the wave function (1) written down 
by Liang et at, to the doped case. We will derive © from a microscopic treatment of the t — J model, based on a 
reformulation of the model in the so-called phase string representation. 

The wave function differs from the Gutzwiller projected BCS wave functions in a fundamental way. Spin 
correlations with various length scales emerge self-consistently as a function of doping in |RVB) and their interplay 
with the charge degrees of freedom is incorporated through the nonlocal phase string operator e~ lB ° in c ia . We will 
find that the ground state @ describes both an AF insulator (at zero doping) and a pseudogap phase with short range 
AF correlations (at finite doping) in which two distinct kinds of spin excitations emerge. We will identify a hidden 
ODLRO in this pseudogap ground state which ensures a finite pairing amplitude and nodal fermionic quasiparticles. 
Although superconducting phase coherence of this state is prevented by the presence of free spins which behave like 
vortices, there exists an inherent superconducting instability at finite doping, which can naturally give rise to <i-wave 
superconductivity. 

The outline of the paper follows. In Section II, we reformulate the t — J model using the phase string representation. 
This is an exact reformulation and allows us to sort out the effects of hole motion in an AF background discussed 
earlier. A singular gauge transformation automatically ensures that the single occupancy constraint is satisfied. A 
brief review of this formalism and important physical consequences will be discussed in Section II. In Section III, we 
obtain the ground state within a mean field theory and in Section IV we transform the ground state wave function back 
to the electron representation, thus completing the derivation of the wave function (J2J . A discussion of the properties 
predicted by J5J follows and finally a superconducting ground state is obtained. The final section is devoted to a 
general discussion about the region of validity (in the phase diagram) of this approach, as well as the questions that 
this study raises. In Appendix A, some detailed derivations of the mean-field theory in the phase string representation 
are presented. In Appendix B, we discuss briefly, an effective theory of spin and charge degrees of freedom that are 
coupled by the phase string. In Appendix C, we show that the ground state obeys time reversal and spin rotational 
symmetries. 



II. PHASE STRING FORMULATION AND PHYSICAL IMPLICATIONS 

In this section, we reformulate the t — J model using the phase string representationii. This representation, like 
the slave-particle representations used in the literatur o 12 ! 13 , is an exact reformulation. However, the phase string 
formulation has a distinct advantage over these other schemes, in the following sense. The most singular effects 
encapsulated in the t — J model are sorted out by a unitary transformation, and all nontrivial phases present in the 
original model now occur as topological gauge fields. Without these gauge fields, the model becomes trivial, i.e., free 
of sign problems, which is the case in one dimension, or in the particular case of half filling in two dimensions. The 
central issue in two dimensions is then how to handle these topological gauge fields at finite doping. It turns out that 
they are well controlled, at least for small doping. 



A. Phase string representation 

The most singular effect of the hole moving in an AF background is the phase string effect. This leads to a 
competition between hopping and superexchange processes. Since it is most singular at length scales of a lattice 
constant (as Z2 signs), it is crucial to track this effect explicitly, before constructing any effective theory. As shown 
by Weng et ai, this can be achieved starting from a slave- fermion representation and by performing a unitary 
transformation of the basis states^. In this procedure, the electron operator is decomposed into 

c la = h\b la e 1 ^ (3) 
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where holon h\ a and spinon bi a operators are both bosonic fields, satisfying the constraint 

>4/n+5>U- = 1 - ( 4 ) 

The phase factor e l ® i(T in Eq. J3J is defined by 

e i&i« = e ie~ (5) 

where TV/, is the (total) holon number operator, the phase string operator Q s i " ng = I [$| _ $0 _ <j<I>'- 1 ] is a nonlocal 
operator with 



(6) 



and 



*? = $>(')> (7) 



^=E^(o^. (8) 



Here, n| Q and nj 1 are spinon and holon number operators respectively, at site I, and &i(l) = Im In [z% — Zi), where 
zi = xi + iyi is a complex coordinate on the lattice. 

It is easily verified that the fermionic statistics of c J(T is automatically ensured by e 1 ®™ " ; besides ensuring fermionic 
statistics, the string operator also incorporates the singular effects of hole motion in the RVB background, as will 
be discussed later. The factor (a) h guarantees anticommutation relations between opposite spins. The factor 
(— a) 1 = ±1 in 10) is added for convenience, and incorporates the Marshall sign into the decomposition (0) of the 
electron in terms of the holon and spinon operators, which may be regarded as a bosonization scheme for the electron 
moving in an RVB background. The same bosonization decomposition can be also obtained from the slave-boson 
representatio n 14 ! 15 . 

B. The t — J model in the phase string representation 

Rewriting the t — J model using the electron decomposition we get Ht-j — H t + Hj, where the hopping term 

H t = -t E ('• , ) h\hj (f" A %) b) a b ia + h.c. , (9) 

and the superexchange term 

"/ tE {^f^i ■ (10) 
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(ij) 



The RVB pair operator Af^ in Hj is defined by 



^ = 5]e^Wj- t ■ (11) 



The unique feature in this formalism is the emergence of three phases (link fields) : , , and , defined on the 
nn links. Without these, there should be no nontrivial sign problem in H t -j, because h and b are both bosonic fields. 
The matrix elements of Ht-j would then be real and negative-definite in the occupation number representation of h 
and b. Consequently, the ground state expanded in terms of these bosonic fields will have real and positive coefficients. 
Therefore, any nontrivial signs (phases) of the t — J model arise solely from these three link fields, defined by 



Afj = - ^ Ml) -9,(1)} [ }_*nl ) , (12) 
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<t>% = \ > 

and 

It is easy to see that the Hamiltonian H t _j is invariant under U(l)xU(l) gauge transformations 

hi -» fce*** , j4|j -» Afj + - Vi ) , 

and 

-> ft^e**' , 4 -> ^ + (0, - ty) . (16) 

Thus andA^ are gauge fields, seen by holons and spinons respectively, as the latter carry their gauge charges 
according to l|15l) and iJTBJ. But they are not conventional U(l) gauge fields, as they must satisfy the topological 

constraints: J2c^-fj ~ n J2iec ( n l1 ~ n ii) ana - J^c^-ij = n J2iec n i according to and ljTl)l . respectively. Here the 
subscript c denotes an arbitrary closed loop such that the fluxes enclosed, ^ c Afj and J2 C ^ij > are determined by the 
number of spinons and holons respectively, inside the loop c; viz., Afj and are not independent gauge fields with 
their own dynamics. Rather they are directly connected to matter fields as a pair of mutual Chern-Simons gauge 
fields^. The term "mutual" refers to the fact that Afj describes quantized ir fluxoids attached to spinons, coupled 
to the holons in H t . Conversely, A^ describes quantized 7r fluxoids bound to holons, coupled to spinons in Hj. The 
phase <j)®j describes a uniform and constant flux with a strength ir per plaquette: ^ D ^ = ±7r according to l|13|) . 

In the case of the one dimensional model, A\j = 0°. = A^ = 0. Thus there is no sign problem in the phase-string 
representation of the t — J model. It implies that the Hamiltonian can be treated within a mean field framework, 
which indeed results in a correct Luttinger-liquid solution for the large-U Hubbard modelii. In two dimensions, these 
nn link phases can no longer be "gauged away" . They compose a mutual gauge structure which completely determines 
the essential sign problem of the t — J model. These gauge fields are generally well controlled in the regimes of our 
interest: 4>% is a non-dynamic phase describing a constant ir flux per plaquette; Afj is canceled when spinons are 
RVB paired at low-temperature phases; A\a remains weak at small doping - it vanishes at half filling, where there is 
no sign problem in the Hamiltonian. Therefore, these gauge fields in the phase string representation are presumably 
well tractable in low doping and low temperature regimes. This is one of the main advantages of the phase string 
representation over other approaches. 

Finally, it is important to note that the phase string representation is defined in a Hilbert space where the total S z 
is an eigen operator. The total numbers of f and j, spinons are conserved individually, such that the topological gauge 
field Afj behaves smoothly as defined in Q12[). depicting fictitious ir fluxoids bound to spinons. (Non-conserved S z 
would result in non-conserved + and — fluxoids in Afj). Different S z states are connected by the spin flip operators, 
defined in the phase-string representation as 

&W- ( 17 ) 

(a factor (— l) Nh has been dropped for simplicity) and S~ — (S^)^ , and Sf = Ylcr a ^ia^i^- These definitions follow 
from (PJ)- The nonlocal phase ^ in l|17f) plays a crucial role in restoring the spin rotational symmetry. 

C. No double occupancy constraint 

Since the t — J model is a projective Hamiltonian, the no double occupancy constraint is a central issue. It imposes 
a severe restriction on the Hilbert space 

Ei^<i, (is) 

which ensures that the system is a Mott insulator at half filling and a doped Mott insulator away from half filling. 
In conventional slave particle decompositions, the constraint reduces to a single occupancy constraint like and 



(13) 



(14) 



(15) 



(-l)V** 
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is still difficult to treat. Typically, the constraint is relaxed at the level of mean field theory, and restoring its effect 
beyond the mean field theory (say, by a gauge theory) is a very challenging problem. In contrast, the phase string 
representation outlined in the previous subsection provides a new way to handle the constraint. 
A wavefunction tp e in the electron c-operator representation is defined in a quantum state as 

|*> = J2 Mi**h ik}) cl lT i T • • • 4 M AAi ■ ■ ■ c k- M Jo> , (19) 

{iu}{j d } 

where the j spin electron sites, {i u } — - ■ -,im, and the J, spin sites, {jd} — ji,- • -i3n c -m, obey the no double 
occupancy constraint in the restricted Hubert space of the t — J Hamiltonian. On the other hand, in the phase string 
representation, the wave function ipb{ii, • ■ ■, U/; ji, ■ ■ -,jN s -M', h, ■ • •, In h ) is defined as 

l*> = E M{iu};{3d};{h})\{iu};{jd};{h}) , (20) 

{iuHid} 

where the indices {lh} = h, • • •, In h denote the empty sites, that are not independent from {i u } and {jd,} under the 
constraint. Here the spinon-holon basis in the phase string representation, {jd}] {lh}} is given by 

\{iuh {j4; ik}) = {-i) N KAi ■ ■ ■ b LAi b li ■ ■ ■ b L-ui h ^ ■ ■ ■ ^.-m 10} , (21) 

where the vacuum |0) is chosen to be filled by holons. The sign factor (-1)^ in can be identified with the 
Marshall sign, and N\ denotes the total number of f spins in sublattice A. Here and henceforth, we will always use 
i to specify an j spin, j a J. spin, and I, a holon, where the subscripts u, d, and h label the sequences of the f spins, 
I spins, and holons, respectively. 

Now, the decomposition (|3J) of the electron operator c^, relates the wave functions written in the two representations: 



where 

z'a . „ zj — Zi , -. r Z7, — z' A 



ud K - z u 1 u \i k - 1 ltd, te* - z 3«> 1 if: \«* - z ^ 1 

= c-'Hizt-zD n«-<^n(4-^) n k-^i fn«-oni*j--^i) ■ 

ud u<_u' dKd' h<h' \ uh dh / 

The coefficient C is given by 

c = n \z lu - z jd 1 yi \ziu- 1 n ~ 1 n ^ ~ ^ \ n ^ ~ z ^ i n _ ^ i 

u<u' d<d' h<h' uh dh 

= J\ \ zk _ Zm i - 

in which k and m run through all lattice sites, such that C is a constant. 

Correspondingly, the wave functions in the electron and phase string representations are related by 

i> e (h, ■ ■ -,«m; ji, • • - Jn c -m) = JC ipb(h, ■ ■ -,«m; ji, ■ ■ - Jn^-m'Ju • • -JnJ , (24) 

which holds generally for the t — J model. 

Note that the above expression for IC has been obtained strictly under the no double occupancy constraint. But 
if the constraint is relaxed in (|23[1 . while still treating C as a constant, then one finds JC — for double occupancies 
when any two particles (spinons and/or holons) occupy the same site; viz., JC defined in (|23|l automatically enforces 
the constraint (J3J through the Jastrow-like factors. Therefore, as far as ip e is concerned, the no double occupancy 
constraint in ipb or ijjb/C is no longer important, since JC in 1)24(1 naturally plays the role of a projection operator. 

Clearly, an exact rpb would satisfy the single occupancy constraint l|4"|. But, the point is that an approximate tpb 
determined without the constraint does not affect tp e) owing to JC. It means that in the phase string representation 
the constraint is indeed "unimportant" , which may be understood in the following way. In the phase string 
representation, the effect of JC in the original tp e , is transformed into the topological gauge fields, Afj and A^, in 
the Hamiltonians, 10 and l|l(J|l . which describe spinons and holons as mutual vortices, as perceived by each other 
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(cf. the discussion in the previous subsection). This clearly implies a mutual repulsion between two species, since a 
spinon cannot stay at the center of its vortex (which is a holon), and vice versa. Thus the constraint that a holon 
and a spinon cannot occupy the same site is now reflected in the interactions present in the new Hamiltonian, and 
the condition Q is not needed as an extra condition to enforce. Note that the constraint (0J also requires the hard 
core conditions among the holons or spinons themselves. But since both holon and spinon fields are bosonic fields, 
local hard core exclusions usually do not involve the sign change of the wave function. Hence, in the phase string 
representation, the local constraint l@J is neither crucial nor singular, as far as low energy physics is concerned. 

D. Phase string effect: A singular doping effect 

In this subsection, we point out that the other important effect, i.e. , the phase string effect, is also explicitly built 
into this representation via IC. 

Let us first rewrite IC (|2*5)l in a more compact form, 



where 



and 



K, = JG , (25) 



j = n ~ z h ) n ^ z u - z *w ) n.( z L - z t ) n k - ^ i n ~ ^ \ n ~ ^ \ > ^ 

u<u d<d' ud h<h' uh dh 



uh |Zl " Zlhl 

It is easily seen that the Jastrow-like factors in J enforce the single occupancy constraint discussed previously: J 
vanishes if two spins (or holes) occupy the same site, or if a hole and a spin occupy the same site. The factor J also 
explicitly captures the fermionic statistics of electrons. 

Now let us focus on the additional factor Q in IC, which is asymmetric with regard to f and j spins: it only involves 
an | spin complex coordinate z* and a holon coordinate z\ h . Suppose that the holon is taken along a closed loop 
shown in Fig. 1. The IC factor will generally acquire a nontrivial phase through Q . In the following we examine this. 

For a spinon inside or outside the loop, there will be no net contribution to Q as the acquired phase is either 2w 
or under such an operation. So any nontrivial contribution is solely from spinons occupying lattice sites on the 
holon's loop-path. At each step of hopping, the holon under consideration will have to switch positions with a spinon 
occupying a nn site on the loop, obeying the no double occupancy constraint. Then Q picks up a minus sign from 
the Jastrow factor in Q every time when the holon and an \ spinon are exchanged. Thus, if one changes the holon's 
coordinate continuously through the loop without encountering other holons, Q will acquire a sequence of signs, 

0— » (-1)^0, (28) 

where A^ op denotes the number of f spinons exchanged with the holon during the above thought experiment. (Note 
that the spinons displaced by holon hopping during the operation should all be restored to the original configuration 
by pure spin flips in terms of the superexchange Hj, and this does not generate any further signs in Q). 

Therefore, the electron wave function %jj e , acquires a sign factor (which can also be thought of as a Berry phase) 
through Q: 



(-iyvp= n(-^)=n^- ( 29 ) 

loop loop 

This factor (for the closed loop) is just the phase string shown in Fig. 1, where a m denotes the index of spins on the 
loop exchanged with a holon, under the thought experiment of moving the holon around a closed path (note that a 
closed loop on a bipartite lattice always involves an even number of nn links). Equation (|29|l shows that the singular 
hopping effect is explicitly incorporated in ip e through /C, in the phase string formulation of the t — J model. 

To conclude this section, the two singular effects in the t — J model, viz., the single occupancy constraint and the 
phase string effect, are explicitly represented by the factor IC, (|23() . Consequently the ground state wave function in 
the phase string representation, -0b, is presumably no longer singular, justifying a mean field theory. 
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III. CONSTRUCTION OF THE GROUND STATE IN PHASE STRING REPRESENTATION 

A mean field theory (bosonic RVB theory) based on the phase string representation of the t — J model was 
constructed in an earlier paper—, and we will follow a similar route here in dealing with the spin degrees of freedom. 
However, an important distinction in this work is the way the charge degrees of freedom are incorporated, in order to 
better facilitate and optimize the kinetic energy of the holes, as well as to construct a state in which the no double 
occupancy constraint is strictly enforced. As shown in Sec. IV, this construction directly results in a pseudogap 
ground state characterized by a finite pairing amplitude and the emergence of nodal quasiparticles. The ground state 
has an inherent superconducting instability which occurs, when phase coherence is realized. 



A. Spin degrees of freedom and bosonic RVB order parameter 

Let us begin by considering the superexchange term Hj, which is expressed in terms of an RVB pair operator Afj 

m (HJ. Note that Af- is invariant under the gauge transformation 1(16(1 . It is then natural to define the bosonic RVB 
order parameter by 

A s = (Af.) (30) 

for nn sites, i and j. This RVB order parameter is different from the RVB order parameter introduced by Anderson 
and collaborators^: it involves the charge degrees of freedom through the gauge field, A^j, and is gauge invariant 

under l|lo|l. Note that \b\bja) = in this theory, and A s is the sole order parameter. At half filling, A\- — in 

\ / Tin 

Af-, and A s reduces to the well known Schwinger-boson mean field order parameter ^2 a (bi a bj— (T )^ [in our definition 
bi a differs from the conventional Schwinger boson representation by a sign factor (—a) 1 ]. The mean field solution 
describes both long and short range AF correlations of the Heisenberg model quite well, as will be discussed below. 
Thus, the bosonic RVB theory we present is ensured to yield good results at half filling, and provides a natural 
generalization away from half filling. 

We now show how the spin degrees of freedom are treated within the bosonic RVB theory. The superexchange 
Hamiltonian is rewritten as 

H. = -J, £ (>^) b\ a b]_ a + h.c. - A£(£ *>lM +5-1), (31) 

(ij)a- i a 

where J s = JA s /2, and A s is assumed to be real and positive. The Lagrangian multiplier A is introduced to implement 
the single occupancy constraint Q on the average. We reemphasize that H s is not a mean field Hamiltonian in the 
conventional sense. Spinons are coupled to a gauge field A^, related to the charge degrees of freedom, and the 
Hamiltonian is invariant under the internal gauge transformation ((16(1 . So H s is still a gauge model, describing low 
energy, long wavelength spin fluctuations, underpinned by the bosonic RVB order parameter 1(50) 1 . The link field A\, 
defined in l|14|l describes holons as 7r fluxoids perceived by spinons in 1 (31) 1 . In the ground state, A^ will be treated as 
a constant field, as though the holons are Bose condensed (which condition will be determined self consistently, when 
we discuss the charge degrees of freedom) , with 



£4-™* (32) 



for each plaquette, where S is the hole concentration. 

Then, H s can be diagonalized by a Bogoliubov transformation 17 



as 



bia = £ W ma {i) (umlmo ~ VmlL-a) - (33) 



H s = £ E m ^ ma 'y ma + const. (34) 



V2 n / x \V2 



Here, u m = f ^ + lj , v m = ^ - lj sgn(£ m ), and E m = ^/X 2 - The Lagrangian multiplier A is 
determined by enforcing ^V (b ia bi a ) = (1 — S)N. The wave function w ma and the spectrum £ m are determined 
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by the following eigen equation, 



£ro™m<r(«) = ~Ja E ^ aA ^ W ma (j) . (35) 
j=nn(i) 



The ground state of H s is constructed in the usual way by imposing the condition, 7 TOC t|RVB)mf = 0, 

|rvb) mf = cxp nr; u ,,/-;.'/ j |o) , (36) 

where the RVB amplitude Wij is given by 



J2—w m i{i)w* ml {j) . (37) 



In deriving the above, we use w^^ii) = w m — a (i), which follows from l|35(l . Furthermore, it is easy to show that for 
each state labeled by to, there is always a state labeled by to, for which = — £ m and w 7 n a (i) = (— l)*u> TOtr (i), such 
that 



1 - (-1) 1- -? 

w« = — V^^' 



7^ only if i, j S different sublattices , (38) 

i.e., the RVB amplitude Wjj only connects j and J, spins on opposite sublattices. 
The nn spin correlations are determined by 

(Si-S J -> nn = -||A'| a <Q (39) 

which is antiferromagnetic in nature [i.e., in the resulting bosonic RVB theory, (Si ■ Sj) > (< 0), if i and j belong to 
the same (opposite) sublattice]. Therefore, A s ^ generally characterizes a regime with short range AF correlations, 
present in a wide range of temperature in the t — J model at low doping (e.g., ~ J/fcs ~ 1, 500 K in the zero-doping 
limit). Fig. 2 sketches a phase diagram in which A s or short-range AF correlations controls the high-energy, short- 
distance correlations in a low-doping regime (< xrvb)^- Such a spin singlet-pairing regime, characterized by A s ^ 0, 
defines the "pseudogap" regime of the bosonic RVB state, with a characteristic temperature T , as illustrated in Fig. 
2. 

In the following two subsections, we further examine the spin correlations at half-filling and finite doping within 
this bosonic RVB description. 



1. Antiferromagnetism at half-filling 

At half filling, the bosonic RVB ground state |RVB) = Pq |RVB)mf, which satisfies the no-double-occupancy 
constraint (Pq denotes the Gutzwiller projection operator), can be explicitly written as 

l RVB > = E (Ell ^ 1 (O*) • ■ • { b Na N h N ) 1°) > (4°) 

{a s } \pair (ij) J 

where each site is occupied only by one spinon, with {a s } = a%, <T2, • • •, &N being partitioned into RVB pairs connected 
by the RVB amplitudes, Wij. Here each RVB pair is denoted by (ij), where i(j) refers to the position of an f (j) 
spin, as before. The summation runs over all possible RVB pairing partitions. The holon annihilation operators in 
(14011 arc introduced since the spinon vacuum |0) is defined as a state where every site is occupied by a holon, for later 
convenience. 

| RVB) given in (|40|l is equivalent to the class of variational wave functions (cf Eq.JTJ) proposed by Liang, Doucot, 
and Anderson (LDA)&. The RVB amplitudes can be determined from the bosonic mean field theory outlined above. 
At half filling, A' 1 - = 0, and the mean field Hamiltonian (|31|l reduces to the Schwinger-boson mean field Hamiltonian. 
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FIG. 2: Schematic phase diagram of the b-RVB mean field theory: The bosonic RVB order parameter A s underpins a pseudogap 
phase at T < To, in which the AF state, spontaneous vortex phase (a low-T pseudogap state), as well as the superconducting 
instability occur at low temperatures, separated by some critical dopings, x c and xbvb, respectively. 

In this case, the RVB amplitude oc |ry |~ 3 , for |ry| ^S> a (a being the lattice constant), as shown in the inset 
of Fig. 3. We use the loop-gas method employed by LDA to determine the nn spin correlation (which also decides 
the ground state energy of the Heisenberg model), as well as the staggered magnetization, m. The maximum sample 
size is 64 x 64. We obtained (Si • Sj) nn = —0.3344(2) J, and a staggered magnetization, m = 0.296(2). These results 
are essentially the same as the best variational result obtained by LDA, and also compare extremely well with exact 
numerical results^, (Si • Sj) nn = —0.3346, and m = 0.31. 

Therefore, at half filling, the insulating state |RVB) correctly reproduces both short range and long range AF 
correlations of the Heisenberg model. AFLRO appears naturally in the long wavelength limit, with a correct magnitude 
of magnetization; short range spin correlations are also well described as evidenced by the excellent agreement of our 
results with the ground state energy obtained from numerical studies. 



At finite doping, the RVB amplitudes Wij, determined by i|3T|) and l|j r 7|) decay exponentially for large spatial 
separations, 



as shown in Fig. 3. Correspondingly, AFLRO disappears and a finite length scale, £, for equal time spin correlations 
emerges; viz., a featureless RVB spin liquid with short range correlations evolves continuously from an AF ordered 
state at half filling (£ — > oo). 

Spin excitations can be created by the action of a spin raising/lowering operator on the spin liquid |RVB) in (|40l) . In 
the bosonic mean field theory, spin excitations are described by the spinon RVB Hamiltonian, (|31|l . The results from 
mean field theory have been reported elsewher e 17 ! 18 , and we briefly summarize the results below. At small doping, 
the spin excitations evolve (from the gapless magnon mode at half filling) into a sharp resonance at low energies, 



2. Spin liquid at finite doping 



|Wy|oce 



(41) 



with a characteristic length scale 




(42) 
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with a characteristic energy scale, E g oc SJ, around the momentum (tt, tt). The spectral weight of this excitation 
is proportional to S, while its width in the momentum space is proportional to which follows from l|41ll . The 
evolution of the gapless magnon mode at half filling into a resonant mode at finite energies reflects the phase string 
effect on the spin degrees of freedom, through the effect of the gauge field, . 



B. Charge degrees of freedom: optimizing the kinetic energy of doped holes 

At finite doping, the mean field state |RVB)mf describes a condensate of RVB spinon pairs, and does not involve 
the charge degrees of freedom (holons) explicitly; i.e., the distribution of holons is independent of the positions of 
spinons. Since the single occupancy constraint is no longer crucial in the bosonic RVB theory, it is tempting to posit 
a mean field ground state 

|GS) = |*>„ ® |RVB) MF (43) 

where \^)h describes the holon condensate. This approach, adopted in a previous worhii, may be justified for a 
long wavelength, low energy theory as the local constraint is not crucial in the phase-string formalism 19 . However, to 
obtain the ground state wavefunction, it is desirable to employ a construction where the single occupancy constraint is 
satisfied exactly; viz., a method that treats spin backflow effects more carefully, which is also important for facilitating 
and optimizing the kinetic energy of doped holes. In the following, we show how this can be done differently. 

First, we project out all double occupancies in |RVB)mf, an d define |RVB) as given in ^4Up. where each site is 
occupied only by one spinon, as if the system were at half filling. However, the RVB amplitudes, Wij, will be 
determined by the bosonic RVB mean field theory of H s at a fixed hole concentration. 

Next, for every doped holon at site I, a spinon at the same site I is annihilated, by the action of h\bi a on |RVB), 
i.e., 

(h\b lrT ) |RVB) (44) 

such that the constraint is automatically satisfied. As discussed below, the construction (|44|l optimizes the kinetic 
energy of the holon-spinon composite, h)bi a , by maximizing the overlap between the RVB configurations before and 
after hopping. Clearly, the behavior of the spinon characterizes the spinon backflow accompanying the motion of 
the holon in the RVB background. 
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FIG. 4: Ratio of the effective hopping integral obtained from 1491 . 

Thus, the full ground state may be constructed as 

\^ G ) = Y,MhM,--^N h )X{Z^h){h\ h b lhah )\mB) , (45) 

{h} h 

where fh{{}h\) is a bosonic (symmetric) many body wave function for holons, and Z ah (lh) denotes a single particle 
wave function for the spinon backflow bi h(7h , with the spin index ah under the constraint l/^X^/i' 7 '* = S z (note that 
S z = in |RVB)). 

In the following, we describe how to determine fh({^h}) and Z ah (lh) in <|45|l based on a mean-field approximation. 
For this purpose, we introduce a mean-field ground state as follows 

|* G ) = Pg\^g)mf (46) 

where the mean-field state 

|*g)mf = ^^(ii,I 2 ,..,k)I]^(y (h\ h bi h a h ) |RVB) MF (47) 
{In} h 

By optimizing the kinetic energy of the holon-spinon composite, h]bi a , based on (|47|) . <fih({lh}) and Z ah (lh) are 
determined as follows. As shown in Appendix A, one finds 

<Ph(h,h,-jN h ) = const. (48) 
namely, the holons are in a Bose condensed state; Z ah (lh) satisfies the following eigen equation 

(-t h )Z ah (() = -i Yl e-^i-^Z^j) , (49) 

j=nn(i) 

with th > being the maximal eigenvalue such that the total hopping energy is given by 

(* G |ff t |* G ) M F/(*G|*G)MF = - tt h SN , (50) 
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and the "bare" hopping integral 



iA*r 

2n b 



n b = 1 - 6 . (51) 



Note that A s ~ 1 at small (5, one has i ~ i in terms of (|51|) . So the holes acquire a good "bare" hopping integral 
in the bosonic RVB background since the singular frustrations due to the phase string effect are gauged away in the 
phase string formalism [one would get a minus sign in the second (RVB) term in (|51|l in the original slave-fermion 
representation] . The residual frustrations of the phase string effect are represented by the phase field 0? • and in 
the eigen equation l|49l) which determines the renormalizcd hopping integral th- The latter is reduced from t as shown 
in Fig. 4 numerically. In Appendix B, we briefly discuss how to go beyond the above mean-field construction in an 
effective model description, which is essentially the same as obtained in Refd 



IV. PSEUDOGAP AND SUPERCONDUCTING GROUND STATES IN ELECTRON 

REPRESENTATION 

So far, we discussed the reformulation of the t — J model using the phase string representation, where the singular 
effects arising from the constraint on occupancy and the phase string, are sorted out explicitly. This allowed us 
to develop an approximate theory in the above section. In this section, we will rewrite our results in terms of the 
underlying electronic degrees of freedom, and discuss the physical consequences and implications. 



A. The ground state (1451) in electron representation 

It is straightforward to reexpress the ground state l|45[l in terms of the electron operators based on © 
|*g) = Y,^{hM,-^N h )J{z ah {l h )(e-^ h c lhah )\RVB) 

{lh} h 

= J2Mli,h,...,lN h )Y[Z ah (l h )£ lhah \RVB) (52) 
{l h } h 

where we introduce a new "hole" creation operator q ct = h\bi a which is related to the electron operator by 

c w = e- i ®*"ci a . (53) 

The "twisted" hole operator, c; CT , satisfies bosonic statistics because of the phase factor, e~ tBi " (see Sec. IIA). On 
using the holon condensation condition (|4*5|) and noting that cf a = and Q|C;JRVB) = 0, the ground state can be 
rewritten in its final compact form 

£z T (/)c ZT J f 5^(06*4 J |RVB) (54) 

where Nj t + N f [ = Nh is the total number of the holes. 

We can further reexpress |RVB) in the electron representation after some straightforward manipulations. By using 
igSJ and ©, we find 

|RVB) = const. ^2 $hvb(oi, ct 2 , ' ' ' 0^)4^4^ ' ' ' c n* n 1°) ( 55 ) 
where the RVB wave function, $rvb, is defined by 

pair (ij) 

Note that |RVB) is always half-filled, involving spins on the whole lattice which satisfy the no double occupancy 
constraint. Here each spin configuration {cr s } = a±, 02, • • •, ctjv is weighted by ^rvb in a way that spins are all paired 
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up with the RVB amplitudes, (— X) l Wij. The latter depend upon doping as determined by the mean field theory 
described in the previous section (see Fig. 3). The summation in ^rvb runs over all possible pairing partitions for 
each configuration {a s }. 

The ground state |^g)> reduces to |RVB) at half-filling, where it describes the antiferromagnetism extremely well 
as discussed in Sec. Ill Al. Upon doping, |RVB) may be regarded as an insulating RVB vacuum. From (|54() . we see 
that each doped hole is "twisted" into a bosonic object c, via the nonlocal phase operator Oi CT , in order to optimize its 
kinetic energy in such an RVB vacuum. The form of the state I'Fg) suggests a generalized "spin-charge separation"; 
viz. , the bosonic holes created by c, move freely in the background |RVB); they carry the charge while the vacuum is 
described by the pairing of neutral spins (spinons). These two degrees of freedom are indeed decoupled in the mean- 
field sense. But this is different from the typical spin-charge separation because the twisted holes are accompanied by 
spins. Indeed, annihilating a spinon a at the site / in the RVB vacuum by ci„ will always create an unpaired spinon 
—a within a scale £ where the spins were originally RVB paired (with no free spinons). So the state \^g) is composed 
both of RVB paired spinons and twisted holes as the holon-spinon pairs as in (57) below. 

Since the spins at the hole sites I's in the vacuum |RVB) are automatically annihilated by Q ff 's, the no double 
occupancy constraint is exactly imposed in \^g) a t finite doping. As shown in Appendix CI, the ground state retains 
time reversal symmetry, which may not be obvious from the form of I'Fg)- R also respects global spin rotational 
symmetry under the approximation that the twisted holes are perfectly Bose condensed, as discussed in Appendix 
C2. 

The structure of \^g) in the doped case looks rather simple in terms of the twisted holes moving on the RVB 
vacuum. But it would look more complicated when written in terms of the c-operators in l)54|l because 0j ff is highly 
nontrivial. Physically, the phase string 0j CT can be regarded as the nonlocal phase shift induced by doped holes. In 
the one dimensional case, where it is responsible for the Luttinger liquid behavior, the phase string can be identified 
precisely with the phase shift obtained from the Bethe ansatz solution for the large-?/ Hubbard modelii. In the 
two dimensional case, the role of the phase string operator 0^ may be seen more clearly in the wavefunction form 
discussed in Sec. II, i.e., (|24[) . Below we rewrite the wavefunction l|54|) in terms of the underlying electrons. 

Starting from the wavefunction ipb, corresponding to the ground state l|45|l in the basis of (|21|l . and using the 
relation between the wavefunctions in the two basis sets Q24[) . we get, 



ipe ({««}; {3d}) 



pair 



(57) 



Here, we have invoked the holon condensation condition l|48|l . In 1)54(1 . a given spin-hole configuration, {i u }, {jd} 
and {lh}, is partitioned into pairs, denoted by (ij) (the sites, i and j, are occupied by an up and a down spin, 
respectively), (hji) (sites l\ and j± are occupied by a hole and a down spin, respectively), (i\l'i) (two sites, i\ and 
l L , are occupied by an up spin and a hole, respectively), and (11') (a pair of holes). The latter three pairings have 
amplitudes, X hn = Z^(lx)Wi xix , X ilVi = (-1)* 1 Wi^Z^h), and Y w = Z 1; (l)Z l (l r )Ww , respectively. Fig. 5 shows 
two possible partitions for a given spin-hole configuration, and the summation in (|57l) runs over all possible partitions. 
Thus, the wavefunction l|57|l shows that the phase string operator 0i CT leads to the generalized Jastrow-type factor K, 
defined in l|25|l . The remanining part of the wavefunction is bosonic, describing RVB pairing of spins, pairing between 
spins and holes as well as between holes. 



B. Hidden ODLRO in the ground state: A pseudogap state 

The properties of the RVB vacuum |RVB) were discussed in Sec. Ill A. In the following we shall focus on the charge 
degrees of freedom in \^g)- 

^,From the definition of l|54|) . it follows that 

(¥ G | (j2 z S)^j I*g) - (58) 

for finite doping: n1" — 0(SN/2). Thus the ground state exhibits an ODLRO as the twisted holes, described by Cj CT 
[or more precisely the bosonic object, Z a (l)ci a , which is gauge-invariant under (|15|1 and (|16l) ]. are Bose condensed 
on the RVB vacuum |RVB) in l|54|) . Since Ci a is a nonlocal operator which differs from the electron operator by an 
infinite-body phase-string twist, the corresponding ODLRO in l|58|l is unlike the conventional ODLRO based on a local 
combination of electron operators, such as spin magnetization or the BCS superconducting pairing operator. In other 
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FIG. 5: Pairing among spins (full circles) and holes (open circles) are illustrated by the connection of lines. Two sets of bonds 
(solid and dashed lines) represent two different resonating pairing patterns, forming a loop covering of the lattice for a given 
spin-hole configuration. 



words, it does not have a direct correspondence to a physical observable. Such a hidden ODLRO characterizes how 
doped charges (holes) are collectively reorganized and behave in the RVB background, |RVB). As the consequences, 
the ground state and excited states will exhibit some very unique novel properties. This will be the subject of our 
discussion below. 



1. Pairing amplitude 

We first show that the hidden ODLRO in \*S?g) Wlu result in a finite pairing amplitude between electrons in the 
singlet channel. 

The electron singlet pair operator can be expressed as 

IT 

— ^ CifjC ~~ i Cj—(j 

= A? J .e i K*l+*J) . (59) 

In the last line, the pair operator is formally written in terms of an amplitude and a phase part, with the amplitude 
operator defined by 20 

A° = £ [(-iye-^-<-^] c^ a (60) 

a 

(a constant factor (— l) Nh is omitted). 
We now demonstrate that (A9-\ ^ 0, i.e., the pairing amplitude is finite as the consequence of the hidden 



ODLRO. Since the twisted holes are Bose condensed, (ci a Cj- a y \^g)n s ^ 2 = A\ 1 $g)n c + excited states, where the 



coefficient A is nonzero and can be estimated by A ~ Z a (i)Z- a {j)\j N^N^ at small doping 2 ^. The pair amplitude 
=n c -2 (* G |A?-|* G )jv e is then given by 



2 

T 1 ! ~ , ATI 



^E^-^^W' ( 61 ) 



where we have used ~ Nfr ~ N h /2 = 8N/2, and 

Z*_ a (j) = Z a (j)(-iye^ . (62) 
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The above equation can be obtained from l|49l) . and by noting that e 1 ** - **^ = — e l2< ^. In order to estimate the size 
of the pair amplitude, we focus on the pairing amplitude between nn sites, (A^) ?m . Invoking translational invariance 
of the ground state, we get 



A = — V (A ) 

<ij> 



<ij> u 
2t 



= j5 . (63) 

Therefore the ground state (|54|l indeed possesses a finite pairing amplitude as a direct consequence of the hidden 
ODLRO in I^g): which is proportional to the doping concentration for small doping. Note that a finite pairing 
amplitude exists only in the singlet channel. It is easily shown that the mean value of the pairing amplitude in the 
triplet channel vanishes identically [in this case, an additional sign a appears inside the summation of 16U|) ]. 

However, superconducting phase coherence is generally absent. The phase part of the pair operator (|59|l involves 
$f which is defined in © and is closely related to the spin degrees of freedom. More precisely, each unpaired spin 
contributes to a free ±2tt vortex via <&f in the pair operator (|59|l . As will be discussed in subsection C, the presence of 

free spins in the ground state \^g) prevents phase coherence in <|59")l . Nevertheless, 7^ in (|5*9T) is a meaningful 

definition and description of a low temperature pseudogap ground state, called the spontaneous vortex phased. In 
this phase, each unpaired neutral spin (spinon) always carries a ±27r supercurrent vortex, known as a spinon-vorteon^. 
We postpone further discussion of the spinon-vortices to Sec. IV C. 



2. Fermionic nodal quasiparticles 

In the following, we construct a quasiparticle excitation based on the ground state I^g)- In particular, we show 
that the existence of nodes in the excitation spectrum is ensured by the hidden ODLRO. 

Consider the action of c< on the ground state \^g)- Since Ci a = Ci a e t&i " and the twisted bosonic holes are 
condensed, one finds 

4l*G>*.-i = e-^°cl\*G)N e -i 

~ e"'2^B„|* G )iv e + excited states , (64) 

where 

B ia ~ v^ + I^We-^e^***) 

^ ^ {a) N h{ _ a y z ^ e i^ e -if(^) . (65) 

Let us focus on Bi„ first. Rewriting Bi„ using l|62l) . we get 



B 



n - v — *l 



SJiziMM-**) , (66) 



and 



% = Vf {-V-ZMe-'M-*!) . (67) 

Nontrivial spatial oscillations in Bi a mainly originate from the phase factor e 2 v ; i > . We rewrite this phase factor 
as a string operator, 

e -**[*?-*f] = e -*f [*?-*f] e i *[*fi-<] e - i *[*?i-<]e i *K-*^] • • • , (68) 
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where i, i±, 12, • •■ is a sequence of nn sites on a path ending at site i. Using the relation — ^(i) — ±ir, and 

invoking the Bose condensation of the holes, we get 

e -*f[*?-*f] e **K-*f 1 ] ; s ie Tfcrfa-^ e -<«'K-<] j etc., (69) 

with 77 = S such that 

B ia Bt a = £>£ e -**o-(r«-r,) (?0) 

where 

ko=(±|[l-^],±|[l-„]) (71) 
at four points in the Brillouin zone (setting the lattice constant a = 1). The quasiparticle weight is proportional to 

~ 0(*) , (72) 

which is gauge invariant, as can be seen from 1)49(1. Here, Fy denotes a shortest-path connecting sites i and j. 
Therefore one finds 

4|* G >tf._i =s e-^^e-'W^N.. (73) 

If phase coherence is realized such that 

e-^*I|* G )~|* G > (74) 

(see Sec. IV C), then cj^ 1^0)^-1 would have a finite overlap with |$G)jV e j namely, the quasiparticle spectrum has 
four nodes at momenta, ko's, that approach (±§, if) in the limit of zero doping. Similar momentum structure can 
be identified in Cko-I^G^c+i- A true d-wave-like quasiparticle state in the superconducting state can be constructed 
as a linear combination of cj_ J^c?) and C-^-a]^ g) > which should be infinitesimally close to the ground state \^g)n c 
at the four nodal points, ko's. 

However, as discussed in Sec. IV C, the phase factor e _I 2*i precludes phase coherence, and therefore the quasi- 
particle state c\ (J \^g)n c - 1 actually has a vanishing overlap with the ground state \^g)n^ because of such a phase 
factor. Thus, the low-lying fermionic quasiparticles near the four nodes, ko's, have a vanishing spectral weight 
Zk = |j\r e ( , I'G|c k(T | v I'G)jv, ! -i| 2 = 0; viz-i the non-superconducting (pseudogap) ground state \^>g) is a non-Fermi-liquid 
in nature. Long lived quasiparticle excitations are expected to emerge only when superconducting phase coherence is 
realized in (|74fl . 

We emphasize that the nodal structure of quasiparticles results from the ODLRO of the bosonic twisted holes, 
in which the phase string factor e~ %Bi ° , reflecting the fermionic nature of the electron operator, determines the 
momentum positions of four nodes. Finally we note that in the previous effective phase string theory in terms of 
holons and spinons (see Appendix B), a similar d-wave quasiparticle excitation can be obtained as the consequence 
of forming the holon-spinon bound pairs by considering residual interactions beyond the effective model^S. 

3. Two types of spin excitations 

The hidden ODLRO in (|54|l is not only responsible for the existence of a finite pairing amplitude as well as d-wave 
like nodal quasiparticle excitations, it also predicts the presence of a new type of low-lying spin excitations which are 
quite distinctive as compared to the previously discussed spin excitations created in the RVB background |RVB) (Sec. 
Ill A). 

The latter may be called Type I spin excitations, which can be created by acting a spin operator, say, Sf , directly 
onto the insulating spin liquid |RVB), without involving the hole condensate in (|54|l . Obviously, this type of spin 
excitations is the only elementary excitation at half-filling, where it reduces to a conventional spin wave in the 
long-wavelength limit as discussed in Sec. Ill A. At small doping, the majority of spins still remain RVB paired, 
characterized by A s 7^ in |RVB). Thus, this kind of spin excitation will naturally persist into the pseudogap regime, 
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where a spin gap will be opened up in its low-energy spectrum as the spin-spin correlations are short-ranged (here Wij 
becomes exponential decay with the distance \i — j \ , see Fig. 3). In fact, as mentioned in Sec. Ill A, the corresponding 
low-lying excitation reduces to a sharp resonance-like mode in the pseudogap regime, with a characteristic energy 
scale Eg ~ 5 J around the momentum (%, n). Remnant spin wave like dispersion is still present at high energies^, 
which disappears when the temperature reaches To as shown in Fig. 2 where the short-range RVB pairing vanishes 
finally. 

Now let us discuss a new type (Type II) of S = 1 excitation which is associated with the twisted holes and may 
be regarded as composed of a pair of quasiparticle and quasihole, instead of a pair of neutral spinons in |RVB). To 
create a Type II spin excitation by, say, S^, a pair of quasiparticles will be created and annihilated in |\1/ G ), which 
involves a spin flip in the "hole condensate" part in (|54H . 



To show this, let us rewrite Sf' = cL.Cjj by 



(-l)V*?] 4^1- (75) 

Then, by using (• ■ -)hole-part *° denote Type II spin excitations created by in the condensed "hole" part in I^g), 
we obtain 



(^ + l*G» holc _ part - (-l)V*- Z t (i)Zl(i)^N^Nl + 1|* G ) . (76) 

The term on the right-hand side (rhs) has the same origin as the quasiparticles around nodal points discussed 
previously. But the phase factor e~ ieia in the quasiparticle channel [ lj64(l ] is now replaced by (— l) l e 4 * i in the 
spin channel, after a recombination of e i<9iT and e~ lQi i in 175|! 24 . The leading contribution to the spin-spin correlation 
at large distance is then given by 

(# G |S7S+|* G ) hole _ part ~ (-I)*-* {z;{j)Z^)Jl%Z^i)Zl{i)) SN (77) 
(using N% ~ y/SN) where 

nfc=e<*f-*?). (78) 

The momentum structure can be determined by the Fourier transformation of (|77|l . in a way similar to that in the 
quasiparticle channel. Note that, an additional topological effect in e 1 ^ (not in e 8 ?*^) requires a little more careful 
treatment. As detailed elsewhere^, (— l) l e l *> h will lead to the following incommensurate momenta, to leading order 
of approximation, 

Q s = [±tt, ±tt(1 - K S )\ and [±tt(1 - K S ), ±n] (79) 

with Ks — 2(5. 

On the other hand, the dynamics of l|77|l is mainly decided by the vortices in e** ; which is in turn determined 
by the charge density. To leading order approximation, no superexchange interaction in the spin degrees of freedom 
influences J7JJ. If one neglects the fluctuation in e I$i , the term on the rhs of (|76|l would have a finite overlap with 
|\I/ G ) at the incommensurate momenta Qa, which implies a zero mode or the presence of an incommensurate AFLRO. 
But the true long-range order will be prevented by the density fluctuations, because any small fluctuation in the 
charge density will be accumulated in the phase factor e I$i in a sufficiently long time and distance, which in turn 
will result in a slow-fluctuating incommensurate AF ordering above I^g); instead of a true static incommensurate 
AFLRO. This low- lying spin fluctuation can provide an alternative explanation, without invoking a dynamic stripe 
picture, of low-energy results in neutron scattering experiments on the cuprate superconductors. 

Therefore, Type II spin excitations in l|76() are quite distinct from Type I in nature. As a direct consequence of the 
hidden ODLRO, they are absent at half-filling [the weight of this term is proportional to the density of the charge 
condensate, from (|76|) ]. One should thus properly distinguish these two in a neutron scattering measurement, as a 
unique prediction of this theory. 



C. Superconducting phase coherence 



As described above, the hidden ODLRO in \^>g) leads to a nonzero pairing amplitude. But it does not necessarily 
imply superconductivity. Consider the action of the superconducting pairing operator on the N e electron ground 
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state. The resulting state has an overlap with N e — 2 electron ground state, determined by using (|59|) and l|63[l . as 
follows: 

( Mf] |* G )^ e ~ A°e- l ^ + ^)\^ G ) Nc _ 2 + excited states. (80) 

V J J nn 

So superconductivity is established only when 

/ e -4(*f+*!)^ , (81) 
viz., the onset of superconductivity is determined by phase coherence. 

1. Spinon vortices 

We noted earlier that each isolated (S — 1/2) spin in |^g) means the creation of a ±2ir vortex in the phase of the 
superconducting order parameter A?p, via according to l]8Up. Thus, each spin automatically carries a phase or 
supercurrent vortex, called a spinon- vortex, in the state described by |^g)- This was first identified in previous work 
based on the effective theorySLS^. 

The ±2-7r vortices centered at f and J, spins in |RVB) are paired as vortex-antivortex pairs, since the underlying 
spins are always paired in the singlet RVB vacuum. Then, one has 

(RVB|e- 4 K*;'+^)|RVB) ^0. (82) 

This is in analogy with the two-dimensional XY model below the Kosterlitz-Thouless transition where vortex- 
antivortex pairs are bound together. In our case, vortices and antivortices are associated with | and | spins, re- 
spectively. The condition that S z — in the singlet background, |RVB), is equivalent to the charge neutrality 
condition in the XY model. Since | and J, spins are RVB pairs with a finite range £ at finite doping [see l|42|l ]. the 
cancellation in ^ + <I>|) is ensured at large distances. Thus, we see that the spin RVB pairing in |RVB) is very 
essential for phase coherence to occur. Note that this condition fails in the AF phase, where £ — * oo, or beyond the 
doping concentration xrvb where the RVB pairing disappears. 

However, the existence of free spins, associated with doped holes created by q ct , will imply the presence of free 
vortices in \^g)- To see this, following (|54|l and (|82|) . one can write 

e-**(^)|* ) = h£ m )c lt J l^Z^cA e-^(W)|RVB) 

„ e^+^c) (83) 

where A^ 1 = J2i^a ^*(0 ~ ^nu w hh denoting the number operator of twisted holes of spin a at site Then, 

the Bose condensation of these free vortices in \^g) will disorder the phase of the superconducting pairing order 
parameter according to l|83|) (although A^ 1 itself vanishes on average, its fluctuations will satisfy an area law in the 
Bose condensed case). 

Therefore the ground state \^g) is not a true superconducting ground state due to the lack of phase coherence, 
although it possesses a pairing amplitude as well as (incoherent) nodal quasiparticles. It is characterized by the 
presence of unpaired spinon- vortices as a direct consequence of the hidden ODLRO in I^g)- Nontrivial Nernst 
effect contributed by such unpaired spinon-vortices has been investigated based on the effective phase string theory 
elsewhere2i. 



2. Superconducting ground state 

The superconducting instability will take place when those unpaired spins associated with twisted holes also become 
RVB paired. It reduces the superexchange energy when those backflow spins form the RVB pairs, and most importantly 
it removes the logarithmic-divergent energy associated with free vortices as vortices-antivortices are paired up at 
low temperature. Consequently the phase is cancelled out in l|83(l and the true phase coherence (|81|l of the 
superconducting order parameter can be finally established. 
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Corresponding to the pairing of those free spins associated with twisted holes, the twisted holes themselves are also 
effectively paired and the ground state should evolve from l|54|l into the following form 



\^g)sc = const. [D 



|RVB) 



(84) 



in which 



^2z a (i)Z- a (j)d 



la^j—a 



(85) 



with g(i — j) characterizing the pairing of the twisted holes. Note that in the limit of no pairing, with g(i — j) = 
constant, equation Q84|l will recover the pseudogap (spontaneous vortex phase) ground state (|54|l at = = Nh/2. 
According to J53J, one has (D)^O. By using (|55|> - l|6"2l . D can be further written as follows 



D ~ J2g(i-j) 



— V, 

6N 



■ —iaA 



A - 

,Zl{i)Z a {j)^- 



(86) 



such that 



(b) = \ E 9(i i)A°. (e-*i(*? + »J) Af) * (87) 

Since the phase coherence ^e - ' 5 ^^*^^) ^ is established self-consistently once g(i — j) introduces the pairing 
among twisted holes as discussed above, one finally obtains (^Afp^ ^ according to Ij87(l . Note that the phase factor 

e _I 2(* i+ *j) also determines the d-wave symmetry of (^Afp^ as discussed previously in Refi^fi. 

Therefore the superconducting ground state l|84ll can be obtained as the result of the pairing instability of twisted 
holes in the pseudogap ground state (|54J) . Simultaneously the phase coherence condition (|81|l is achieved, and fermionic 
electron operators in D become meaningful as 

D\^ G )sc ^Y.9i3Z~2 aCl " c 3-^ G ^ sc ^ ( 88 ) 

ij a 

with the d-wave symmetry characterized by 




Equation (|88|l also implies the restoration of the coherent fermionic quasiparticle excitations in the superconducting 
state. 



V. SUMMARY AND DISCUSSION 



In this paper, we obtained and analyzed the ground state of a pseudogap phase (or spontaneous vortex phase) 
for a doped Mott insulator, starting from a generalized mean field treatment of the t — J model in the phase string 
representation. In this representation, the most singular effects of doping, viz., the no double occupancy constraint, 
and the phase string effect arising from the disordering of the Marshall sign, are sorted out through a systematic 
procedure. All nontrivial phases in the phase string representation are characterized by a topological gauge structure, 
which turns out to be well controlled in the regime of interest (low doping). The physical picture of the ground 
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state that emerges, is that of an underlying RVB vacuum (a spin liquid insulator) which makes the doped holes be 
twisted into bosonic objects, thereby gaining kinetic energy. The ground state retains time reversal and spin rotational 
symmetries, and has a simple form (|54[l in the electron second-quantization representation. 

At half filling, the ground state describes an AF Mott insulator, and provides a very accurate description of both 
long and short range spin correlations of the Heisenberg model. At finite doping, the spin degrees of freedom are 
described by an RVB spin liquid l(55|l : the mean RVB pair size (f^2"|) reduces continuously as doping increases. Charge 
degrees of freedom emerge at finite doping, with a hidden ODLRO l|58|) . This leads to a finite (electron) pairing 
amplitude, which characterizes the pseudogap structure of the ground state. The hidden ODLRO in the charge sector 
is also found to be the underlying reason for the four nodes in the quasiparticle excitation spectrum, but the ground 
state is non-Fermi-liquid like in the absence of phase coherence. It is more appropriate to call the pseudogap ground 
state as a spontaneous vortex phase, since free S — 1/2 spins that act like vortices prevent superconducting phase 
coherence in this state. 

Two types of low-lying spin excitations above the pseudogap ground state were identified. One is composed of 
neutral spinons and has a resonance-like feature around (tt, tt). It reduces to the conventional gapless spin wave 
at half-filling. The other type of spin excitation exists only in the doped regime and describes a slowly fluctuating 
incommensurate AF ordering. It provides an alternative explanation of the neutron scattering experiments without 
invoking a charge stripe scenario. 

Superconducting phase coherence is established by an instability which causes RVB pairing of the free spins asso- 
ciated with the twisted holes. The ensuing ground state (|84|l . is non-BCS-like in many aspects. The properties of the 
superconducting ground state (|84|l . especially those involving non-BCS-like features, are very interesting problems to 
pursue further. 

The class of wavefunctions discussed in this paper provides a unified description of Mott physics, antiferromagnetism, 
pseudogap physics, and d-wave superconductivity within a single framework. There are many basic distinctions 
between this ground state and BCS-like wave functions. An important distinction is the independence of the spin 
excitations from the quasiparticles. The presence of the pairing amplitude and spontaneous vortices associated with 
spin excitations is another unique feature and is independent of superconducting phase coherence. Furthermore, 
generalized Jastrow-like factors K. in (|23|l that appear in the ground state wave function, are incompatible with the 
(Slater) determinantal wave functions of the BCS theory. Therefore, the pseudogap and superconducting ground 
states obtained in this work represent and describe a new state of matter. 

Given that the state 1|54|) undergoes a superconducting instability with phase coherence at low temperature, we 
may discuss further implications for the global phase diagram. Clearly, the phase coherence condition l|82|) fails in the 
spin ordered phase, when the spatial separation of the RVB pair, £ — > oo; indeed, at half filling, the RVB amplitude 
Wij shows a power law decay at \i — j\ — > oo, without a finite length scale, as discussed in Sec. Ill A. In this AF 
phase, phase coherence in (|82|) due to the tight binding of vortices and antivortices in <I>f is disrupted. In general, 
such a spin ordered phase will persist into finite doping at sufficiently low hole concentration, S < x c 26 . The condition 
of phase coherence is not satisfied also when the RVB pairing between spins in |RVB) is continuously reduced to 
zero, at sufficiently large doping xrvb- In this case, the underlying rigidity for binding of vortices and antivortices 
in the ground state is also destroyed and the superconductivity should be suppressed beyond xrvb- How the ground 
state wavefunction systematically evolves with the doping concentration and the possible existence of quantum critical 
points at, say, x c and xrvb, is an important subject for further investigation. 

Finally, the results presented in this paper can be used as a starting point to pursue several interesting questions, and 
we indicate a few possible directions here: (i) The pseudogap (spontaneous vortex phase) ground state is inherently 
unstable to superconductivity as the result of pairing of the twisted holes. Nonetheless, the pseudogap ground state 
can still be regarded as a true ground state, for instance, when strong magnetic fields are applied to destroy the RVB 
pairing of the low-lying spins associated with the backflow of holes. The study of this state in the presence of strong 
magnetic fields is thus experimentally relevant, (ii) It will be very useful to develop a numerical scheme based on 
the wavefunction (|oTjl for a variational study. At half-filling, we used the wavefunction to compute the ground-state 
energy, magnetization, and equal-time spin-spin correlations by employing the loop gas method. As pointed out, 
this gave us accurate results. However, away from half filling, the RVB amplitude Wij is no longer a real function 
and loop products generate a fictitious flux proportional to doping concentration [cf. <|C28|) ]. This leads to a sign 
problem and further work is necessary to overcome this barrier. Havilio and Auerbach studied correlations in doped 
antiferromagnets using a combination of numerical methods and the Gutzwiller approximation for a different class of 
RVB wave functions^. A similar study of the wave function proposed in this paper should yield interesting results, 
(iii) Finally, we note that the form of the wave function l|24|l follows from very general consequences of a hole moving 
in an RVB background, since the disordering of the Marshall sign, i.e., the phase string effect, occurs in systems 
even without AF long range order. Thus the phase string representation and the bosonic RVB mean field theory can 
be used to study models such as the doped quantum dimer model, where the spin background is known to be RVB 
paired. 
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APPENDIX A: DETERMINATION OF ip h AND Z a (l) 
/From the mean- field state 11471 . one finds 



;mf 



(*g|*g)mf = J2 l^({^})| 2 (RVB|nz;,(/, l )Z ffh (^)(v<) (bl h< bi h a h ) |RVB) r 

!i<i2<— h 

« E \<Ph({k})\ 2 ]l\Z, h (l h )\ 2 {RyB\nl a jRVB) MF . (Al) 

ii</2<--- h 

In obtaining the last line, we omit RVB pairings involving two holon sites, e.g., the second term on the rhs of the 
following expression (note that lh refers to a holon site): 

(RVB« t <JRVB)mf = (RVB \n\^\ RVB) M f (RVB |nf 2 J RVB) M f + (RVB|&J iT &^RVB) M F(RVB|fe ilT & i2i |RVB) M F. 

The essence of this approximation is that for a dilute concentration of holes, the hole- hole correlation induced by the 
background spin RVB pairing can be neglected, since the average separation between holes sets the upper bound for 
the size of the RVB pair wave function. For the sake of clarity, we shall invoke this approximation in determining 
(Ph({lh}) and Za- h (lh) below. In principle, these effects can be incorporated without affecting our results qualitatively. 
Now, consider the hopping term Q, which is rewritten as 

Ht = -tJ2 (e^-^-^) (hK) (hjbl) + h.c. (A3) 

It can be interpreted as describing the hopping of a holon-spinon composite h\bi a under the influence of the gauge 
field A*j — — crAfj in the RVB background |RVB); b ia represents the spinon backflow accompanying the hopping 

of the holon h\. 

The | and J, spinons are paired in the RVB background |RVB), and it costs finite energy to break up an RVB pair 

at finite doping. Correspondingly, (Af^) = and (^A\^) 2 ^ remains small in (|A3|) because Afj depicts fictitious ±tt 

fluxoids bound to a — ±1 spins. We point out that this condition could fail both, at very low doping, (6 < x^), where 
a spin ordered phase exists, and at large doping (5 > zfwb), where the RVB pairing disappears and A\^ fluctuates 
strongly. In both limits, we expect different phases. 

Neglecting A\^ and the hard core repulsion among holons, one may simply choose 

(ph(h,h,.» Jnh) = const. (A4) 

as if the holons are in a Bose condensed state. Consistent with this picture of holon condensation, A\a can indeed be 
treated as describing a constant flux ttS per plaquette. We recall that such an approximation [cf. I|32|l ] was made in 
determining the RVB amplitudes, Wij. 

We can now determine Z afe (i^) by minimizing (>£ c\Ht\^ g)mf / '( v t 'g\^ g)uf- Upon invoking the approximations 
detailed above, we get 

* <; '-'EE ^-^ZkW^j) , (A5) 



<*g|*g)mf h=1<ij> 



where 
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By optimizing the hopping energy under the condition \Z ah (i)\ = 1, we get the hopping energy 

(* G |tf f |* G )MF/(*G|*G>MF = ~ tthSN , (A7) 

with th > being the maximal eigenvalue of the eigen equation 

(-t h )Z„ h {i)=- t l Y, e-^~^Z ah U) , (A8) 
j=nn(i) 

and Z„ h , its eigen wave function. 

^From (| A6|) we see that t has two contributions: the first term, proportional to (1 — S)/2, is the probability that 
the hole hops between nn sites, in a background of uncorrelated spins; the second term, proportional to |A S | 2 , is 
the enhancement from RVB assisted hopping - the initial state (before hole hopping) and the final state have a non 
vanishing overlap owing to the RVB pairing in the spin background. Note that t ~ t at small 8 as n b ~ 1, A s ~ 1. The 
effective hopping integral th can be obtained from i|A8(l . for a fixed hole concentration and under the approximation 
(|32ll . The numerical results are shown in Fig. 4. 

Finally, we note that the eigen values of (|A8|) would form a Hofstadter spectrum, if the gauge field A\^ were treated 
at a mean field level, as a uniform flux l|32|) . However, the density of the backflow spinons is tied to that of the holons 
and so the fluctuations in due to fluctuations in the holon density are self consistently related to the density 
fluctuations of the spinon backflow. Thus, (|A8|) resembles an anyon (semion) system in the bosonic representation 
with a spin index 27 . The mean field approximation is usually no longer valid in determining its excitation spectrum. 
This is a nontrivial issue. Nevertheless, since we shall be only interested in its ground state, which is a condensate 
with a uniform density distribution, all we need to know is that there is a renormalized constant th for the effective 
holon hopping term, estimated as the maximal eigen-value at the mean- field level in (£EJ. The wave packet Z a (i) is 
determined as a linear combination of the wave packets of a cyclotron radius a c = ci/VttS in accord with (|32[) and is 
made maximally uniform in space, in order to accommodate the condensation of the holes. 

APPENDIX B: EFFECTIVE THEORY 

In the bosonic RVB mean field theory, the holon many-body wave function tph has been assumed to satisfy the 
ideal Bose condensation condition (j48(l at T — 0. As we noted in the paper, two effects are omitted in l)48[l. viz., 
the hard core correlation among holons and the fluctuations of the link field Afj . The latter effect should become 
important when excited spinons are present. Furthermore, ifh will also become nontrivial in the presence of an external 
electromagnetic field. 

A general form of the wave function ifh in l|47|) under the influence of Afj, and in the presence of an external 
electromagnetic field A\j, can be determined by the following effective hopping term 

H h = -t h J2 e iA ^ +iA *i h\hj + h.c. (BI) 

(ij) 

If the hard core condition of the holon field is neglected, and if A\a = 0, 1|B1|) leads to the Bose condensation solution 
(I48|l in the ground state. Note that Hh is a gauge model consistent with the gauge invariance l|15|l in the original 
Hamiltonian. Hh also respects the spin rotational symmetry as one can easily check that [Hh, S] = 0, where S is the 
total spin operator. Furthermore, in the absence of Af -, the time reversal symmetry of Hh can also be shown, by 
noting that Afj — ► — Af ■ under the flip of the spins. 

The effective hopping integral th appears in (|Blf) as a renormalized t, which is decided by the spinon backflow 
according to (|49|1 . Here the spinons do not directly see the gauge field Afj and the external electromagnetic field j4|,-, 
because they satisfy a different gauge transformation (|16[1 in l|49|l . Thus, the spinons are truly charge neutral and only 
carry S — 1/2 in the phase string formalism, in contrast to, say, the slave boson gauge theories where both holons 
and spinons are coupled to the same external electromagnetic field through the Ioffe-Larkin rule. 

An effective Hamiltonian H stT i ng based on Hh and H s can be written down as 

-^string = Hh + H s . (B2) 

Such an effective Hamiltonian has been derived earliest, and is the basic low energy effective model for the doping 
and temperature regimes underpinned by the bosonic RVB order parameter A s ^ in the phase diagram of Fig. 2. 

We remark that -Hstring in (|B2|I is obtained under the assumption that the spinon backflow accompanying the 
hopping of the holons only provides a renormalized hopping integral th- Here th is determined at the mean field level, 
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as the minimal eigen value of 1|49|) . However, in principle, the internal excitations of the spinon backflow can effectively 
reduce according to 1491) and thus increase the kinetic energy of -ff s trmg- This should be taken into account when 
the high energy part of the charge degrees of freedom are studied. Finally, we point out that the holon and its spinon 
backflow is actually bound together, as previously shownSi by including the residual interactions of the t — J model 
beyond the effective Hamiltonian. Such a twisted hole as a bound holon-spinon pair emerges naturally in the present 
wavefunction approach given in the main text. 

APPENDIX C: TIME REVERSAL AND SPIN ROTATIONAL SYMMETRIES OF THE PSEUDOGAP 

GROUND STATE (??) 

1. Time reversal symmetry 

The time reversal relation for an S = 1/2 single-particle wavefunction is 



Correspondingly, in second quantized language, the time reversal of the electron operators read 



It is easy to check that 

TS.T- 1 = -Si. 
Then the time reversal of the RVB state |RVB) is given by 

T|RVB) = const. $^ VB {a x , a 2 , ■ ■ -a N ) [a x a 2 ■ ■ ■ a N ] 4_ CT1 4_ CT2 ■ • ■ c ] N _ aN |0) 

= const. ^iW- "!;-^,- • ;-^N)cl cri 4 <r2 ■ • -cjv^l ) 
= |RVB) 

by noting that a\a% ■ ■ ■ ajsr = in a bipartite lattice, and 

*rvb(-°'i>- "2 



pair (ij) 



Enc- 1 )^- 1 )^ 

pair (ij) 

EIE- 1 )^ 

pair (ij) 

^RVB(ci, CT 2 , • • -, CTjv). 



(CI) 

(C2) 
(C3) 

(C4) 



(C5) 



(C6) 



Furthermore, by using 



according to Eq.JSJ) and £i a — e i<r CuT, one finds that 
Finally, the time reversal of the ground state |^g) is given by 



(C7) 
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T\*g} = T\^Z^l)~ Cl A [Y^Z^l'YcvA T- l T|RVB) 

= ^^(Ocue-^C-l)*^-!)^ V-^^(Oc n e-^(-l)^(-l) r ] |RVB) 



iV' / v Nt 



= i-i)»k-i) NhlNh - iy2 \£zM«i) (5>(Oq't) IRVB) 

= |*G> 

at iV^ = = Nh/2 (or S z = 0). In obtaining the third line, the condition (|H2"|l is used. 

Thus, we have proved the time reversal symmetry of the pseudogap ground state l|54|). for S z = 0. 

2. Spin rotational symmetry 

The condition for a state with SU (2) spin rotational symmetry is 

S 2 1*> = (C8) 

or equivalently, 

5 2 |*) = 0, (C9) 

5*1*) = 0. (CIO) 

Here S is the total spin and S — S x ± iS v , where S x , S y and S z are the three components of S. In the following we 
shall consider |* G ) with S z = 0, and show that 5 I± |* G ) = 0. 
By noting that, for i ^ I, 

S+e-^t = e -«**fs+ e «iW (Cll) 
S+e<**f = e^'S+ (C12) 



we have 



and 



S+c, a = S+e-* e '°c la 

= e ie ^c l<7 S+(l-6 a ) (C13) 

S + Q ff = Q.5>+e iei(l) . (C14) 



Then according to and , 



Oh} h 

= EII^W^ (^S'+e i S^ e ".«|RVB)] . (C15) 

In terms of (|55|) and l|56[) . one has 

|RVB) = £ *RVb(<ti,<7 2 ,- • -^)cL lC L 2 • • • 4^10) 

= '•'"^•XXU ,: 1: ' u '<./ • (C16) 

{a s } pair (ij) 
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where each spin configuration {oi, 02, ■ ■ -cr/v} is partitioned in pairs denoted by (ij), with i and j belonging to different 
sublattices connected by (— VfWij in i|(J16|) . For each pair of (ij) and (ji), the rest of |RVB) remains the same (the 
nature of RVB pairing) and thus needs not to be considered. Then we find 



Hj 



= ((-ljV^^Wji + (-l^^^Wy) |- • ;i T,i ?,•••) 

= e i (e^m^+e^-mj)) fo. ie ^& - u,, ( ' r ) i- • r,j ?,•••) (ci7) 

where in obtaining the last line, &i h (i) — 9i h (j) — @i(lh) ~ 9j(hi) has been used. 

In the ground state, the holes are uniformly distributed in space due to the ODLRO, which ensures (|32|) . We shall 
show (see below) that under the condition (|32|) one has the following relation, 

Wij = \W ij \e iA % (C18) 
up to a pure gauge transformation. Since | Wij \ = \ Wji | , then one has 

kj = (C19) 

for any given set of RVB pair (ij) in l|C16|) . 

Note that one should be careful about the case that a pair of sites (ij) may be occupied by holes. Obviously there 
is no contribution in (|C15I) if two sites are both occupied by holes, or by a hole and an up spin. If only one of them 
is occupied by a hole and the other is by a down spin, then two configurations, with the hole at site ij j and a down 
spin at site j/i, will have the same amplitude but opposite signs similar to the case in <|C17|) . Since the twisted holes 
are Bose condensed, smearing out the charge distribution one still finds a cancellation with vanishing contribution to 
(|C15ll . Therefore, generally one has 

S+\V G )=0. (C20) 

Proof of fC18\) under the condensation condition I Sty) . According to the definition i|37|l . Wij is determined by l|35|l . 
which is rewritten as 



-,-iA 



j=NN{i) 



e 



"«»mC?) = (C21) 



with w m = w m ^(i) and £' m = —£ m /Js, with a gauge choice 

4<+* = °' ( C22 ) 

= W)** ( C23 ) 

for (|32l) under the "holon" condensation condition. Since the system is translational invariant in the y-direction, we 
may express w m (i) as w m (i) = e lnlyly g m (i x ) — e lmyly g mx , my (i x ), with the eigenfunction (|C21H becomes 



g m (ix - 1) + 9m(ix + 1) + 2 cos(n5i x - m y )g m (i x ) = £' m g m (ix), (C24) 

where m y = 2irn/L, n = 0,1, L — 1. If the doping concentration 8 satisfies ttSL = 2kir (k S Z), the periodi- 
cal boundary condition can be chosen (the Hofstadter case). The more general case can always be infinitesimally 
approached from the Hofstadter case. 

In terms of <|C24|) . we have the following identities 



gm x ,m y (^x ~t~ 1) 9-m x ,m y —7r5 (i x )e lx ™ (C25) 

and 

^m x ,rn y ^>m x .m y — 7r<5' (C26) 

where % m depends on the choice of the phase of g m . Using the above identities, we obtain 

Qm x jUiy^x Tj x ) — Q77i x ,7n y — 7r5r] x (^x)^-^~' Tl— 1 ^ rnJr ^ n 1 ) x m (C27) 
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Substitute the above solution into Wy , 

Wij = Wm {i)w* m {j) 



Then 



Urn im v (iy-j v ) i- \ * (■ \ 

m m 

W i+V>j+V = -^^e im ^- j ^g mx>my (i x +r] x )g* mximy (j x + ri^ 

m 

6 y y y 9m x ,m y -7TSr) x (' l x)9 mxjrny - 7r Sr] n: Ux) 

m Um 

E'^Jimy+irSrixXiy-jy) (■ \ * (a ) 

u m 

m 

and it is straightforward to show 

loop loop 

on a loop (i —> j —> j + T} —> i + 7} —> i). One can always choose a proper gauge such that (|C18|1 holds on any link 
(i-.j)- 
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